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Abstract 

The form factor bootstrap approach is used to compute the exact contributions in 
the large distance expansion of the correlation function < a{x)a{U) > of the two- 
dimensional Ising model in a magnetic field at T = Tc. The matrix elements of 
the magnetization operator a(x) present a rich analytic structure induced by the 
(multi) scattering processes of the eight massive particles of the model. The spectral 
representation series has a fast rate of convergence and perfectly agrees with the 
numerical determination of the correlation function. 


1 Introduction 


Over the past few years, considerable progress has been made in the nse of conformal 
invariance methods and scattering theory for the nnderstanding of the critical points and 
the nearby scaling region of two-dimensional statistical models (see, for instance 0, ^). 
At the critical points, the correlation fnnctions of the statistical models fall into a scale- 
invariant regime and their compntation may be achieved by solving the linear differential 
eqnations obtained by the representation theory of the inhnite-dimensional conformal 
symmetry pj. The sitnation is different away from criticality. The scaling region may be 
described in terms of the relevant deformations of the hxed point actions. These deforma¬ 
tions destroy the long-range flnctnations of the critical point and the associated qnantnm 
held theories are usually massive. If an inhnite number of conserved charges survive the 
deformation of the critical point action, the corresponding QFT can be efficiently char¬ 
acterized (on-shell) by the relativistic scattering processes of the massive excitations. In 
this case, the integrals of motion severely restrict the bound state structure of the theory 
and force the S'-matrix to be elastic and factorizable [|, i- Once the exact S'-matrix 
of a model is known, one can proceed further and investigate the off-shell behaviour of 
the theory by means of the Form Factor approach B 0- This consists in computing 
matrix elements of the local helds on the set of asymptotic states and reconstructing 
their correlation functions in terms of the spectral density representation. As argued in 
[§], and also con&med by the explicit solution of several two-dimensional massive QFT 
[|, 0 0, 0 0, a general property of the spectral series of massive integrable QFT is 
their very fast rate of convergence for all distance scales. This important quality of the 
spectral representation series may be regarded as the key point to the success of the Form 
Factor approach. 

The aim of this paper is to apply the Form Factor approach and compute the spin- 
spin correlation function G{x) =< (j(x)cr(O) > of the Ising model in a magnetic held h at 
T = Tc (in the sequel, this model will be referred to as IMMF). For the other integrable 
deformation of the critical Ising model, i.e. thermal, the spin-spin correlation function has 
been exactly determined in 0,0 and this result, together with its remarkable connection 
with the Painleve’ functions |^, 0, may be regarded as one of the main accomplishments 


in statistical mechanics. On the other hand, very little is known about the spin-spin 
correlation function for h ^ 0 at T = Tc whose determination has been a long-standing 
problem of statistical mechanics. As we will show in this paper, the computation of G{x) 
can now be approached analytically, thanks to the scattering formulation of the model 
proposed by Zamolodchikov [Q. Apart from the importance of computing G{x) itself, 
there are other related theoretical issues which render this calculation instructive. 

The hrst issue concerns the rich structure of bound states and higher-order poles of 
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the S'-matrix. Higher-order poles of the S'-matrix of two-dimensional held theories can be 
naturally interpreted as singularities associated to multi-scattering processes |T^, |T^ ^ 


On this basis, one would expect a similar hierarchy of higher-order poles in the form 
factors as well. There is, however, an important difference between the pole structures of 
the S'-matrix and those of the form factors. In fact, the S'-matrix contains simultaneously 
information about the s-channel as well as the u-channel. Correspondingly, the poles 
of the S'-matrix are always arranged in pairs, with their positions located in crossing 
symmetrical points. This is in concordance with the two different ways of looking at 
a scattering process, i.e. in the direct or in the crossed channel. On the contrary, the 
M-variable plays no role in the calculation of the form factors as these only depend on 
the s-variables of all subchannels of the asymptotic states on which the matrix elements 
are computed. The absence of u-channel singularities in the form factors implies that 
their analytic structure may be different from that of the S'-matrix and it is therefore 
an interesting problem to understand how the higher-order poles enter the form factor 
calculation. 

The second theoretical issue that emerges from the computation of the spin-spin cor¬ 
relation function in the IMMF is a careful reconsideration of the so-called minimality 
prescription, which is usually invoked for computing the form factors. Briefly stated, this 
consists in assuming that the form factors have the minimal analytic structure, compati¬ 
ble with the nature of the operator and the bound state pattern of the scattering theory. 
So far, this prescription has been successfully applied to solve several two-dimensional 
QFT as, for instance, those of refs. gig, 0113,0, despite the fact that the theoretical 
reason of its validity was lacking. One of the basic motivations why minimality is often 
adopted is because frequently the asymptotic behaviour of the matrix elements is not 
easy to determine. In this paper, we will present a simple argument which will allow us 
to place a quite restrictive upper bound on the high energy limit of form factors. Using 
this criterion, one can see that the minimality prescription is violated in the IMMF and 
therefore extra polynomials in the variable s have to be included in the matrix elements 
of the held a{x). These polynomials, nevertheless, can be uniquely determined since the 
scattering theory always provides a sufficient number of constraints. As a matter of fact, 
the equations which £x the extra polynomials are usually overdetermined and this gives 
rise to self-consistency conditions which are indeed fulhlled by the IMMF form factors. 

The last point we would like to mention is the successful comparison of our theoretical 
determination of < a{x)a{0) > with numerical simulations. These simulations have been 


carried out in the last few years by two different groups |2^ and in particular, 

a collection of high-precision numerical estimates of the spin-spin correlation function of 
the IMMF, for different values of the magnetic held and diherent sizes of the lattice, can 
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be found in [^. Although there is no doubt about the actual existence of higher-mass 
particles, which can be easily extracted by transfer matrix diagonalization of the IMMF 
^ or directly from the lattice |]26l, it was quite difficult to see the presence of these 


higher-mass states from the analysis of the numerical data of the spin-spin correlation 
function. Namely, a best fit of the two-point function seemed to be compatible with 
the only contribution of the fundamental particle [^, a result which appears intriguing. 
In fact, from a theoretical point of view there is no reason why decoupling of higher- 
mass particles should occur in the spin-spin correlator since the IMMF has apparently 
no selection rules related to any symmetry. Indeed, as we will see, the magnetization 
field a(x) couples to all particles of the theory and it is only the small values of the 
relative couplings which could be responsible for a possible misleading interpretation of 
the numerical results. 

The paper is organized as follows. In Section 2, we briefly outline the main features 
of scattering theory of the IMMF, the exact mass spectrum of the model and some other 
quantities which will prove useful in the sequel of the paper. In Section 3, we analyse 
the Form Factor approach and the properties of the spectral representation series of the 
correlation functions. In Section 4, a general constraint on the asymptotic behaviour of 
the form factors is introduced and applied to the computation of the matrix elements of 
the magnetization operator (j{x) of the IMMF. We also discuss the occurrence and the 
interpretation of higher order poles in the form factors. Comparison of our theoretical 
results versus numerical simulations, as well as the saturation of the sum-rules satisfied 
by the spin-spin correlation function, are presented in Section 5. The paper also contains 
three appendices. In Appendix A we discuss the general bootstrap approach to the 
computation of the form factors of the model. Appendix B gathers useful mathematical 
formulas to deal with the monodromy properties of the matrix elements. Appendix C 
presents specific examples of matrix elements with singularities associated to higher order 
poles in the scattering amplitudes. 


2 Scattering theory 

In the continuum limit, the IMMF may be regarded as a perturbed CFT. At the critical 
point, (T = Tc and h = 0), the Ising model is described by the conformal minimal model 
Ads ,4 with central charge ^=11^. There are three conformal families in the model, those 
of the identity, magnetization and energy operators, denoted respectively as [/ = (0,0)], 


a = 


and 


e = 


(i, l) , where the numbers (Aj, Aj) in the round brackets are 
their conformal weights. We can move the system away from criticality by modifying the 


(- - 
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action of fhe critical point as 

A = AQ + h J d^xa{x) . (2.1) 

For small values of the magnetic field, the system is still at T = Tc. However, the 
coupling to the magnetic held induces a mass scale M{h) in the problem and destroys the 
long-range huctuations of the critical point. Using the action (|2.1| ), or the corresponding 
lattice Hamiltonian, one can in principle analyse the properties of the model in terms of 
a perturbative expansion in h There are, however, some important questions 

which cannot be easily addressed by using such perturbative formulation, as, for instance, 
the determination of the exact mass spectrum of the model or the computation of its 
correlation functions. In order to answer such questions, we have to rely on a non- 
perturbative approach that exploits the most important feature of the dynamics of the 
model, namely its integrability. The IMMF has, in fact, an inhnite number of conserved 
charges of spin s = 1, 7,11,13,17,19, 23, 29 (mod 30). As shown in a remarkable paper 
by Zamolodchikov the existence of these conserved charges allows one to dehne a 
self-consistent scattering theory which provides an alternative route to the perturbative 
formulation. Let us briefly outline the main results discussed in 0. 

The scattering theory which describes the scaling limit of the IMMF contains eight 
different species of self-conjugated particles Aa, a = 1,..., 8 with masses 


Till 

m2 

ms 

m4 

mg 

mg 

my 

mg 


= Mih), 


IT 


2mi cos — = (1.6180339887..) mi, 

5 

2micos— = (1.9890437907..) mi, 

1 30 V J 1, 

2m2 cos — = (2.4048671724..) mi, 

30 

2m2Cos— = (2.9562952015..) mi, 

15 

2m2 cos — = (3.2183404585..) mi, 

30 

4m2 cos — cos — = (3.8911568233..) mi, 
5 30 

TT 27r 

4m2 cos — cos — = (4.7833861168..) mi, 
5 15 


( 2 . 2 ) 


Within the standard CFT normalization of the magnetization operator, hxed by the 
equation 

^ 0 (2.3) 


< a(a;)cT(0) >- , , 


Fr 



the overall mass scale M{h) has been exactly determined in p3 

M{h) = Ch^ 


(2.4) 
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where 


4 

5 


4sm|r(i) / 4^^r(|)r^(a) \ 

r (!) r (^) V r 0 ) (a) J 


= 4.40490858... 


(2.5) 


The scattering processes in which the eight particles are involved are completely elastic 
(the hnal state contains exactly the same particles as the initial one with unchanged 
momenta) and, due to the factorization of multiparticle scattering amplitudes induced by 
integrability, they are entirely characterised by the two-particle amplitudes Sab (Figure 
1). These are functions of the relativistic invariant Mandelstam variable s = (pa-hPb)^ or, 
equivalently, of m = (pa — Pb)^- Sab has a branch root singularity in the variable s at the 
threshold s = {rria + mb)^. By crossing symmetry, an analogous branch point also appears 
at the threshold of the w-channel, namely at s = {iria — mbY (Figure 2). Those are the only 
branch cuts of the S'-matrix, due to the elastic nature of the scattering processes. The 
other possible singularities of the scattering amplitudes Sab are simple and higher-order 
poles in the interval (m^ — m^)^ < s < {nia + rribY which are related to the bound state 
structure. 

An important simplihcation in the analysis of the analytic structure of the S'-matrix 
comes from the parameterization of the external momenta in terms of the rapidity variable 
6, i.e. = rria cosh 9a, p\ = ma sinh 6a- The mapping s{6ab) = + ml -|- 2mamb cosh 9ab, 

where 6ab = 9a — 9b, or equivalently u{9ab) = siin — Oab) = -|- ml — 2mamb cosh. 9ab, 

transforms the amplitudes Sab into meromorphic functions Sab{9ab), which satisfy the 
equations 


Sab{9)Sab{-9) = 1 ; (2.6) 

Sab{9) = Sab{ni-9) , (2.7) 


expressing the unitarity and the crossing symmetry of the theory, respectively. The simple 
poles of Sab{9) with positive residues are related to bound state propagation in the s- 
channel, as shown in Figure 3, whereas those with negative residues are associated to 
bound states in the w-channel. Suppose that the particle Ac with mass squared 


ml = ml + ml + 2mamb cos Uab, G (0, vr) 


( 2 . 8 ) 


is a stable bound state in the s-channel of the particles Aa and Ab- In the vicinity of the 
resonance angle 9 = the amplitude becomes 


Sab{9 ~ iul^) ~ 


z(F 


ab) 


^ - iKb 


(2.9) 


with F((^ denoting the three-particle coupling^. Since the bootstrap principle gives the 

^Note that the ^'-matrix cannot determine the three-particle couplings but only their square. This 
results in an ambiguity for the sign of all which can be solved by defining a consistent set of form 
factors. 
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possibility to consider the bound states on the same footing as the asymptotic states, the 
amplitudes Sab are related each other by the functional equations Bi 


Sa{e) = S,,{e + iu%) - tui), (2.10) 

= tt — For the IMMF, the solution of eqs. (|^ ), ( p.Tf ) and (|2.10|) are given by 


s'ai,(»)= n (/«(9)r 

O'G.A.ab 


where 


fa{e) ^ 


tanh ^ {0 + ina) 
tanh — ma) 


( 2 . 11 ) 


( 2 . 12 ) 


The sets of numbers Aab and their multiplicities Pa, specifying the amplitudes (|2.11|) , can 
be found in Table 10. In order to correctly interpret this collection of data, notice that 
the functions /a( 6 *) satisfy the equation fa{d) = fi_a{6) as well as /a( 6 *) = faiiT^ — 6). 
Hence, they have two poles located at the crossing symmetrical positions 9 = ina and 
6 = i7r{l — a). Therefore, the poles of the S'-matrix of the IMMF will always appear in 
pairs. 

Applying eq. (p^) to the simple poles with positive residues at 6 ^ = ^, 6 ^ = ^ and 
0 = ^ in S'ii( 6 ') (related to the bound states Ai, A 2 and ^3 respectively), we can extract 


pi — 
111 — 


pa _ 

111 — 


p3 _ 
1 11 — 


\ 


2 tan ^ tan tan 

_ 6 _ 15 


tan ^ tan ^ 

15 10 


^ = 10.990883. 




2 tan ^ tan tan 


15 


30 


tan ^ tan ^ 

15 D 


= 14.322681.. 


(2.13) 


\ 


2 tan ^ tan ^ tan ^ 

--^2-20 ^ 1.0401363. 


tan ? tan 


Stt 


6 10 

Other three-particle couplings can be obtained similarly. In addition to simple poles, the 
S'-matrix of the IMMF presents higher-order poles due to multi-scattering processes^ The 
odd order poles correspond to bound state poles while those of even order do not. Their 
appearance is an unavoidable consequence of the iterative application of the bootstrap 
equations (p.l0|) . In relation with the calculation of the Form Factors, some of these poles 
will be considered in Section 4 and in Appendix C. 


^Note that the numbers Aab of Table 1 should be read in units of 

^ A detailed analysis of the nature of higher-order poles in the S'-matrix of the (1 -t 1) integrable theories 
can be found in [pM ^ . 
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3 Correlation Functions and Form Factors 


Once the spectrum and the S'-matrix are known, we can investigate the off-shell behaviour 
of the theory. We can compute the two-point (as well as higher-point) correlation functions 
of the model through the unitarity sum (Figure 4) 


oo „ 

< $(x)d)(0) >= E / ^ 

n =0 47r 


<O|$(O)|kl,,(0i)---kl„J0„)> 


dOn 

2g-|a:| cosh6»j, 


(3.1) 


Basic quantities of this approach are the Form Factors (FF), i.e. the matrix elements of 
the local operators <F on the asymptotic states (Figure 5), dehned as 

F* .. (3.2) 


A detailed discussion of the form factor approach and the mathematical properties of the 
matrix elements (|3.2|) can be found in |^, |^. Here we simply consider the basic equations 
we need for our discussion. 

For a scalar operator 4>(x), relativistic invariance requires that its form factors depend 
only on the rapidity differences di — Oj. The elasticity of the scattering processes, together 
with the crossing symmetry and the completeness relation of the asymptotical states, 
permit to derive the following monodromy equations satished by the FF 


, (^ 1 ) dii .., 0 jj) 


di ,..,di,di^i , ..dfi ' 


(t' 2 > ■ ■ ■ , t'n, ^l) ■ 

(3.3) 


a2,a3,...,an,ai 


In terms of the s-variable of the channel (Aq. Aaj), the hrst equation in (|3.3|) implies that 
in the form factors there is a branch cut in the s-plane extending from s = {rriai Arrta^Y 
s = oo. This is similar to what happens in the S'-matrix. However, a difference between 
the analytic structure of the FF and the S-matrix comes from the second equation, which 
on the contrary shows that the FF do not have a w-channel cut extending from s = —oo 
to s = {rua^ 

Apart from these monodromy properties, the FF are expected to have poles induced by 
the singularities of the S-matrix. A particular role is played by the simple poles. Among 
them, we can select two different classes which admit a natural particle interpretation. 
The hrst type of simple poles are the so-called kinematical poles related to the annihilation 
processes of a pair of particle and anti-particle states. These singularities are located at 
da — da = in and for the corresponding residue we have 


-thmie-d)Fy. 

e^e 


a,ai, 


^ {d + in, d, 9i 
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FZ .„„(9i,...,9„). 

(3.4) 


This equation can be graphically interpreted as an interference process due to the two 
different kinematical pictures drawn in Fig. 6. 

A second class of simple poles is related to the presence of bound states appearing 
as simple poles in the S'-matrix. 9 = and F^^ are the resonance angle and the 
three-particle coupling of the fusion Aa x Aj, ^ Ac respectively, then FF involving the 
particles Aa and Ai, will also have a pole at 6^ = with the residue given by (Fig. 7) 

-i hm^ .= (3.5) 

where 9c = {9autc + 

In general, the FF may also present simple poles which do not fall into the two classes 
above. In addition, they may also have higher-order poles and indeed, their analytic 
structure may be quite complicated. We will came back to this point in Section 4 where 
the specihc example of the IMMF will be discussed. 

Although eqs. ( |3.4|) and ( |3.5| ) do not exhaust all pole information, nevertheless they 
induce a recursive structure in the space of the FF, which may be useful for their deter- 
minationQ Finding the general solution of eqs. (|3.3|) , (|3.4|) and (|3.5|) for the IMMF poses 
a mathematical problem of formidable complexity, as described in Appendix A. 

Fortunately, because of an important property of the spectral series (0 , an accurate 
knowledge of the correlation functions can be reached with limited mathematical effort. 
This property consists of a very fast rate of convergence for all distance scales |^. In 
view of this, the correlation functions can be determined with remarkable accuracy by 
truncating the series to the hrst few terms only. This statement appears to be obviously 
true in the infrared region (large Mr), where more terms included into the series only add 
exponentially small contributions to the final result. The fast rate of convergence in the 
crossover and ultraviolet regions seems less obvious. In fact, for small values of the scaling 
variable Mr, the correlators usually present power-law singularities and all numbers of 
particles are in principle supposed to significantly contribute to the sum. However, this 
does not seem to be the case of integrable QFT with sufficiently mild singularities in 
the ultraviolet region for a “threshold suppression effect” discussed in |^. Although this 
result was originally derived for QFT with only one species of particles in the spectrum, 
we expect that it also applies to the IMMF|, due to the smooth singularity of G{x) at 
the origin, i.e. G{x) ~ x~^l^. 

^ An important general aspect of the Form Factor approach which is worth mentioning is that the 


validity of eqs. (3.3), ( |3.4D and (3.5) do not rely on the choice of any specific local operator <i)(a:). This 
observation, originally presented in |]^, may be used to classify the operator content of the massive 
integrable QFT, as explicitly shown in refs. ^ ^ 

®As we will show in Sect. 5, this will be indeed confirmed by: (a) a direct comparison with the 
numerical determination of the correlation function; (b) the saturation of the sum-rule related to the 
second derivative of the free-energy of the model (the zero-moment of the correlation function) and (c) 






If this crucial property of the spectral series is taken for granted, the first terms of the 
series are expected to saturate the values of the correlation function with a high degree 
of precision and we can only concentrate on their analytic determination. The number of 
terms to be included in the series essentially depends on the accuracy we would like to 
reach in the ultraviolet region and, to this aim, it is convenient to order them according to 
the energy of the particle states. For the IMMF, the hrst seventeen states are collected in 
Table 2. The most important contributions to the sum come from the one-particle states 
Ai, A 2 and ^ 3 , and for the correlation function we have correspondingly 

I I ^ I I ^ I I ^ 

G{r) = \< 0 |(t( 0)|0 >1^ H- —Ko{mir) H- —Ko{m2r) H- —Ko{m3r) + O 

(3.6) 

where Kq{x) is the modihed Bessel function and we have defined 

T, =<0|a(0)|kl, > . (3.7) 


These matrix elements will be exactly determined in next section. Concerning the vacuum 
expectation value of o'(a;), it can be easily obtained from the relationship between this 
field and the trace Q{x) of the stress-energy tensor. Since a{x) plays the role of the 
perturbing field in the theory under consideration, it is related to Q{x) as 

©(x) = 2ti h{2 - 2A^) a{x) (3.8) 


On the other hand, the vacuum expectation of Q{x) can be exactly determined by the 
Thermodynamic Bethe Ansatz and its value is given by |]3^ 


where 


, ^ / M Trmf 

< O|0(O)|O >=-i , 

Vii 


/ 271 271 71 \ 

(pii = 2 ^ sin vra = 2 ( sin-h sin-h sin — 1 . 

Qe.4ii V 3 5 15/ 


(3.9) 


(3.10) 


Hence, using the above formulas and eq. (|2.4|), we have 

4(^2 

<0|a(0)|0>= -= (1.07496..) 


15(p 
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(3.11) 


Eq. ( p. 6 |) provides the hrst terms of the large-distance expansion of the correlation 
function < (j(a;)cr(0) >. A more rehned determination of G{x) may be obtained by 
computing the FF of the higher-mass states, as discussed in the next section and in 
Appendix A. 

the saturation of the sum-rule derived from the c-theorem (second moment of the correlation function). 
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4 Form Factors of the IMMF 


In the framework of the Form Factor bootstrap approach to integrable theories, the two- 
particle FF play a particularly important role, both from a theoretical and from a practical 
point of view. From a theoretical point of view, they provide the initial conditions which 
are needed for solving the recursive equations. Moreover, they also encode all the basic 
properties that the matrix elements with higher number of particles inherit by factoriza¬ 
tion, namely the asymptotic behaviour and the analytic structure. In other words, once 
the two-particle FF of the considered operator have been given, the determination of all 
other matrix elements is simply reduced to solve a well-dehned mathematical problem. 
From a practical point of view, the truncation of the spectral series at the two-particle level 
usually provides a very accurate approximation of the correlation function, which goes 
even further than the crossover region. This section is mainly devoted to the discussion 
of the basic features of the two-particle FF in the IMMF. 

In the general case, the FF must be a meromorphic function of the rapidity 

difference dehned in the strip Im6 G (0,27r). Its monodromy properties are dictated by 
the general equations (^7^ , once specialized to the case n = 2 


= Sat{e)F!,i-e ), (4.1) 

F^^izTT + 6) = F!,{z7r - 6) . (4.2) 

Thus, denoting by a solution of eqs. (|4.1|) , (|4.2|) free of poles and zeros in the strip 

and also requiring asymptotic power boundness in momenta, we conclude that 
must be equal to F^^^{6) times a rational function of cosh^. The poles of this extra 
function are determined by the singularity structure of the scattering amplitude Sab{d)- 
A simple pole in F^fj{9) associated to the diagram of Fig. 8 corresponds to a positive 
residue simple pole in Sab{9) (see eq. m) and in this case we can write 


Fl{d ^ i<fe) ^ 


,-pc 

ab 

■ - ^ c 


lU 


ab 


(4.3) 


The single particle FF F^ is a constant because of Lorentz invariance. Other poles induced 
by higher order singularities in the scattering amplitudes will be considered later in this 
section. The kinematical poles discussed in Section 3 do not appear at the two-particle 
level if the operator 4>(x) is local with respect to the helds which create the particles, 
which is the case of interest for us. 

Summarizing, the two-particle FF can be parameterized as 

Fl{0) = . (4.4) 
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where Dab{0) and Qabi^^) ci-re polynomials in cosh 0 : the former is fixed by the singularity 
structure of Sab{G) while the latter carries the whole information about the operator <f)(a:). 

An upper bound on the asymptotic behaviour of FF and then on the order of the 
polynomial Q%{0) in eq. (^4.4|) comes from the following argument. Let 2A$ be the scaling 
dimension of the scalar operator <F(x) in the ultraviolet limit, i.e. 


< <h(a;)<F( 0 ) > 


\x 


4A<, 


\X\ 


0 


(4.5) 


Then in a massive theory 


Mp = J (fx \x\P < <F(a;)$(0) >, 


< +CX) 


if 


p > 4A$ — 2 , 


(4.6) 


where the subscript c denotes the connected correlator. The two-point correlator may be 
expressed in terms of its euclidean Lehmann representation as 


< d)(x)$(0) >c= J cPp 


-I- p? 


(4.7) 


where the spectral function p is given by 


1 oo 

^ 27r ^^Xi>...>6i„ 27r 27r 


rlf) dO ^ ^ 

^ • •. ,6'„)|^(5(^mfeCosh6'fc-/i)d(^mfcSinh6'fc) 

k=l k=l 


Substituting eq. (^4.71) into the dehnition of Mp and performing the integrations over p, /i, 
and X, one hnds 


M, ~ Y. 


9\'> ...'>6n 


dOi... dOr, 


. 




(4.8) 


\k=l 


Eq. ( ^4.6| ) can now be used to derive an upper bound for the real quantity y^, defined by 


lim 

0i|—»CXD 


F 


-J- 

ai,...,an 


{9i, . . .,9n) 


rsj 




(4.9) 


This can be achieved by hrstly noting that taking the limit 9i —>■ -|-cx) in the integrand 
of eq. (|4.8|) , the delta-function forces some other rapidity 9j to diverge to minus inhnity 
as —9i. Since the matrix element ■ ■ ■ ,9n) depends on the rapidity differences, 

it will contribute to the integrand a factor in the limit \9i\ —>■ oo. Then eq. ( |4.6|) 

leads to the constraint 

< A$ . (4.10) 

Note that this conclusion may not hold for non-unitary theories since not all the terms in 
the expansion over intermediate states are guaranteed to be positive in these cases. 
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Let us see how the aforementioned considerations apply to the specihc case of the 
IMMF. An appropriate solution of eqs. (|4.1|) and ( |4.2|) , corresponding to the scattering 
amplitudes reported in Table 1, can be written as 


ut”(9) = (-"siniT]" n (Gcwr . 

\ aeAab 


where 


Ga{0) = exp < 2 


dt cosh (a-fjt _ ^{in- e)t 


'o t cosh I sinh t 


sm 


27r 


For large values of the rapidity 

G„(^) ^ exp(|^|/2) , 1^1 ^ oo , 


(4.11) 


(4.12) 


(4.13) 


independent of the index a. Other properties of this function are discussed in Appendix 
B. 

An analysis of the two-particle FF singularities which will be described later on in this 


section, suggests that the pole terms appearing in the general parameterization eq. (|4.4|) 
could be written as 

Dab{0)= n (n-aWr , (4.14) 

where 

ia = n + 1 , ja = n , if = 2U 1 ; 

ia = n , ja = n , if = 2n , 

and we have introduced the notation 


(4.15) 


V4ff) = 


COS ira — cosh 6 
2 cos2 ^ 


(4.16) 


Both F’™”(6') and Dab{d) have been normalized to 1 in 0 = in. 

Finally, let us turn our attention to the determination of the polynomials for 

the specihc operator we are interested in, namely the magnetization held a{x). In view 
of the relation (p.8|) , this is the same as considering the analogous problem for the trace 
of the energy-momentum tensor Q{x). For reasons which will become immediately clear, 
we will consider the latter operator in the remainder of this section. 

The conservation equation = 0 implies the following relations among the FF of 

the diherent components of the energy-momentum tensor 




p ai,. 

..,a„(^li • • 

■■,0n) ; 

(4.17) 



^ —F^ 
P+ “i.- 


■ . , 6'n) , 

(4.18) 
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where ± x^ are the light-cone coordinates and The requirement 

that all the components of the energy-momentum tensor must exhibit the same singularity 
structure, leads to conclude that the FF of 0(x) must contain a factor P+P“. However, 
the case n = 2 is special because, if the two particles have equal masses, the mismatch 


of the singularities disappears in eqs. ([4.18|) and no factorisation takes place. From this 
analysis, we conclude that for our model we can write 


Ql{0) = fcosh0 + "P„,(0), 

\ 2mamb J 


where 


Nab 


p.>.m = i: i 


cosh^ 9 


(4.19) 


(4.20) 


A;=0 


The degree Nab of the polynomials Pab{G) can be severely constrained by using eqs. ( |4.9|) 
and ([4.10|) . Additional conditions for these polynomials are provided by the normalization 
of the operator 0 (x), that for the diagonal elements P^, reads 




(4.21) 


Using all the information above, we can now proceed in the computation of the IMMF 
form factors, starting from the simplest two-particle FF of the model, namely P®( 6 ') and 

First of all, by using eqs. (|4.1CI|) and ([4.13|) one concludes that A^n < 1 and Nu < 1. 
In view of the normalization condition ( [4.21| ), only one unknown parameter, say is 
necessary in order to have the complete expression of P®(0). On the contrary, we need two 
parameters, a °2 and 0^2 > to specify P®( 6 *). To determine all three unknown parameters, 
note that the scattering amplitude S'ii( 6 *) possesses three positive residue poles at 9 = 

9 = and 9 = which correspond to the particles Ai, A 2 and A 3 respectively; on the 
other hand, Si2{9) exhibits four positive residue poles at 6 ^ = 9 = 9 = and 

9 = associated to Ai, A 2 , A 3 and A 4 . Hence, since three poles are common to both 
amplitudes and no multiple poles appear in both of them, eq. ([4.3|) provides a system of 
three linear equations which uniquely determine the coefficients 0^2 and 0^2 


■pc 

11 




= 


pc 

■*- 12 


Resfl=,;„c_ (0) 


c= 1,2,3 . 


(4.22) 


The result of this calculation can be expressed in terms of the mass ratios rrij = rrii/mi 
as 


^TTTTl^ 

Pn{9) = ——(2 cosh6*-I-2-t- 77137717 ) 
77137717 

-^ 12 ( 6 ') = ^^12 (^27772 cosh 6 * -f- 77 I 2 -F 777g^ 


(4.23) 

(4.24) 
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where 


Hi2 = (1.912618..) . 

Equations (|4.3|) can now be used to obtain the one-particle form factors F® (a = 1,..., 4), 
whose numerical values are reported in Table 3. In particular, 


F® = 


vrm 


2 (l-Fmsmr) 


G2G2G± 

3 5 15 V / 


cos^ f cos^ ^ 

5 30 


F® = 


F® = 


r{im3m7 

Airml sin | (2 cos ^ + 2 -f rhsmT^ 


sin ^ sin ^ sin ^ sin ^ 

15 15 10 30 


G2G2G^ , , , 

3 6 15 V 3 /. 


( 27rA 


COS g COS 


i cosf)‘ 


271 


sm ^ 

dvrm^ sin (2 cos ^ + 2 + rh^niT^ 


sin ^ sin sin ^ sin 

5 15 D 30 


G2G2G± 

3 5 15 


(s)] (' 


cos g COS g COS' 


—V 

30 y 


sin ^ sin ^ sin f sin 

10 30 D 30 


rf^rngmy sin 

In order to continue in the bootstrap procedure and compute the other one-particle 
and two-particle FF, we have to hrstly consider the multiple poles of the scattering am¬ 
plitudes. Such poles are known to represent the two-dimensional analog of anomalous 
thresholds associated to multi-scattering processes |^. These are processes in which the 


two ingoing particles decay into their “constituents”, which interact and then recombine 
to give a two-particle hnal state. In the general framework of relativistic scattering theory, 
the location of this kind of singularities is determined by the so-called Landau rules 10 . 
In the two-dimensional case, such rules admit the following simple formulation: singular¬ 
ities occur only for those values of the momenta for which a space-time diagram of the 
process can be drawn as a geometrical hgure with all (internal and external) particles 
on mass-shell and energy-momentum conservation at the vertices. The simplihed two- 
dimensional kinematics only selects discrete values of the external momenta for which 
such a construction is possible and this is the reason why in two dimensions the “anoma¬ 
lous” singularities appear as poles rather than branch cuts. The order of the pole and its 


residue can be determined using the Cutkosky rule pT|] which states that the discontinuity 
across the singularity associated to the abovementioned diagram is obtained evaluating 
it as if it were a Feynman graph but by inserting the complete scattering amplitudes at 
the interaction points and by replacing the internal propagators with mass-shell delta- 
functions 0(p^) 6 {p‘^ — m?). For a diagram containing P propagators and L loops, P — 2L 
delta-functions survive the L two-dimensional integrations; since the singularity whose 
discontinuity is a single delta-function is a simple pole, the graph under consideration 
leads to a pole of order F — 2L in the amplitude . 

Let us initially consider the second order poles. A second order pole at 9 = ip occurs 
in the amplitude Sab{9) if one of the two diagrams in Figures 9.a and 9.b can be actually 
drawn, namely if 

V= - W^cd - Ude ^ [0, tt) . (4.25) 
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The quantity ir] is the rapidity difference between the intermediate propagating particles 
Ac and Ac- From these hgures, it is easy to see that 


V = <d + Udb-'^ ■ 


(4.26) 


The crossing symmetry expressed by eq. (p^.7| ) obviously implies that, in addition to the 
double pole m. 6 = iip, an analogous pole must be present in 9 = i{7i — ip). Since the 
residues of the two poles are now equal, it is impossible to distinguish between a direct 
and a crossed channel, and the two poles must be treated on exactly the same footing. 
At the diagrammatic level, this fact is reflected by the possibility to hnd a diagram 
satisfying eq. ( |4.25|) also for 6^ = z(7r — <^ 9 ). Hence, let us consider only one of these poles, 
the one located ai 6 = iip. In the vicinity of this pole, the scattering amplitude can be 
approximated as (see Fig. 9.a) 


SabiO) ~ 


iKA?ScM 

{6 — iipY 


(4.27) 


Note that the expression of this residue, which is obtained for rj > 0, is also valid in the 
limiting situation rj = 0 (Fig. 9.b), for which a residue (F“^ F^gF^jFjg) is expected. In fact, 
the consistency of the theory requires 


■pa p6 _ p6 pa 

^ cd^ de ^ cf^ fe ’ 


(4.28) 


an equation that is indeed satished for the three-point couplings of the IMMF. Moreover, 
in the case f] = 0, the “fermionic” nature of the particles, expressed by the relations 


SabiO) 


— 1 if a = b ] 
1 if a^b , 


(4.29) 


implies that the two particles Ac and Ag propagating with the same momentum in Fig. 9.b 
cannot be of the same species. In this case, the factor Scdv = 0) in eq.( |4.27| ) equals unity. 

The double pole at 9 = iip in Sab{9) induces a singularity at the same position in 
F^Ij{9). For rj > 0, this is associated to the diagram on the left hand side of Fig. 10. Since 
the singularity is now determined by a single triangular loop, the form factor Flie) has 
only a simple pole rather than a double pole. The residue is given by 


r:drlSce{v)F!c{-v) 


(4.30) 


Eq. ([4.1|) can now be used to write (see the right hand side of Fig. 10) 


Fabi^ ~ i(p) ~ i 


.. ^AFctiv) 


9 


lip 


As written, this result also holds for 7^ = 0. 


(4.31) 
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The poles of order p > 2 in the scattering amplitudes and the corresponding singu¬ 
larities in the two-particle FF can be treated as a “composition” of the cases p = 1 and 
p = 2. This is the case, for instance, of a third order pole with positive residue a.t 6 = iip 
in Sab{0)- In the S'-matrix, a third order pole occurs if the scattering angle r] in Fig. 9.a 
coincides with the resonance angle The corresponding diagram is drawn in Figure 11 
and in this case we have 


Sab{0 ~ i(p) ~ i 




(4.32) 


{6 — iifY 

The third-order pole at the crossing-symmetric position 6 = {{tt — ip) has negative residue 
since it corresponds to the crossed channel pole. With respect to the case p = 2, the pole 
a.t 9 = iip in becomes double (see Figure 12) 


f1{0 ~ ip) ~ 


pa pb p/ 
cd de ec 

{9 - ip)^ ^ 


(4.33) 


while the pole ai 9 = i{7i — p) stays simple. 

The above analysis suggests the validity of the following general pattern for the struc¬ 
ture of the form-factor poles: a pole of order 2n at 9 = ip in the crossing-symmetric scat¬ 
tering amplitude Sab{9) will induce a pole of order n both at 9 = ip and at 9 = i{7i — p) 
in the two-particle form factor FYi,{9)', viceversa, a positive residue pole of order ( 2 n -f 1) 
at 9 = ip in Sab{9) will induce a pole of order {n + 1) at 9 = ip and a pole of order n 
at the crossing symmetrical position 9 = iiji — p) in FYb{9). We have used this result to 
write the parameterization of eq. ([4.14| ). 

Moreover, in integrable QFT, these arguments can be easily extended to the higher 
matrix elements ^^{9i,..., 9n) with n > 2 , since the complete factorisation of mul¬ 
tiparticle processes prevents the generation of new singularities. In other words, the 
singularity structure of the n-particle FF is completely determined by the product of the 
poles present into each two-particle sub-channel. 

We have used eqs. (b:.3|), ( |4.31|) , ([4.33|) to continue the bootstrap procedure for the 
two-particle FF of Q{x) in the IMMF up to the level ^ 3 ^ 3 , and an illustration of the 
method through a specihc example may be found in Appendix C. The results obtained 
for the coefficients (the only unknown quantities in the parameterization of eq. ( |4.4|) 
after the pole structure has been hxed) are summarized for convenience in Table 4; Table 3 
contains the complete list of the one-particle matrix elements. Two important comments 
are in order here. The hrst is that in all the determinations of F^ except F® and F^, 
a number of equations larger than the number of unknown parameters is obtained. The 
fact that these overdetermined systems of equations always admit a solution^ provides a 

^Solutions can only be found by choosing the three-point couplings either all positive or all 
negative. Hence, this restricts the ambiguity of the three-point couplings to an overall ± sign only. We 
are not aware of any other explanation for this constraint on the F^j,. 


16 












highly nontrivial check of the results of this section. The second point is that, since the 
pole structure has been identihed, there is no obstacle, in principle, to continue further 
the bootstrap procedure and to achieve any desired precision in the determination of the 
correlation function in the ultraviolet region. Actually, we will show in the next section 
that the information contained in Tables 3 and 4 are more than enough for practical 
purposes. Nevertheless, from a purely theoretical point of view, it would be obviously 
desirable to have a complete solution of the recursive equations. A possible approach to 
this non trivial mathematical problem is suggested in Appendix A. 


5 Comparison with Numerical Simulations 


The data collected in Tables 3 and 4, together with the vacuum expectation value eq. 
and the three-particle matrix element 62 , 63 ) given in Appendix A provide us 

with the complete large-distance expansion of the correlator < 0(a:)0(O) > up to order 
g-(m2+m3)|x|_ gf |;Pg clggree of convergence of the series, and then of its 

practical utility, is obtained by exploiting the exact knowledge of the second and zeroth 
moments of the correlation function we are considering. Indeed, in a massive theory the 
c-theorem sum rule provides the relation 

3 


C = — y (fx\x\‘^ < 0(a;)0(O) >c , 


(5.1) 


where C is the central charge of the conformal theory describing the ultraviolet hxed 
point. For the Ising model C = ^. In addition, if we write the singular part of the free 
energy per unit volume as fs — —UM‘^{h), a double differentiation with respect to h leads 
to the identity 

U = ^ [ cPx < 0(a;)0(O) >c . (5.2) 

TT^ J 

On the other hand, the exact value of the universal amplitude U is obtained by plugging 


eq. (|3^) into 


U = 


dvr 


< O|0|O >= 0.0617286.. 


(5.3) 


The contributions to the sum rules ( b.l|) and (|5.2|) from the hrst eight states in the spectral 
representation of the connected correlator are listed in the Tables 5 and 6, together with 
their partial sums. The numerical data are remarkably close to their theoretical values. 
Notice that a very fast saturation is also observed in the case of the zeroth moment, 
despite the absence of any suppression of the ultraviolet singularity. 

Let us now directly compare the theoretical prediction of the connected correlation 
function Gc{x) =< (j(x)cr(O) >c with its numerical evaluation. A collection of high- 
precision numerical estimates of Gc{x), for different values of the magnetic held h and 
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different size L of the lattices, can be found in the reference |^. We have decided to 
consider the set of data relative to L = 64 and h = 0.075, where the numerical values 
of G{x) are known on 32 lattice space (Table 7). Such a choice was dictated purely by 
the requirement to use data where the effects of numerical errors are presumed to be 
minimized. Errors can be in fact induced either from the hnite size L of the sample or 
from the residue fluctuations of the critical point, which may be not sufficiently suppressed 
for small values of the magnetic held h. 

In order to compare the numerical data with our theoretical determination, we only 
need to hx two quantities. The hrst consists in extracting the relationship between the 
inverse correlation length, expressed in lattice units, and the mass scale M{h) entering 
the form factor expansion. The second quantity we need is the relative normalization 
of the operator a{x) dehned on the lattice, denoted by aiat(a^), with the operator a{x) 
entering our theoretical calculation in the continuum limit. Let us consider the two issues 
separately. 


The correlation length ^ is easily extracted by using eq. (|3.6|) to analyse the exponential 
decay of the numerical data collected in Table 7. As a best ht of this quantity, we obtain 


M{h = 0.075) = r = 5.4(3) . 


(5.4) 


Let us turn our attention to the second problem. The easiest way to set the normal¬ 
ization of a{x) with respect to cTiat(x) is to compare their vacuum expectation value. The 


lattice determination of this quantity can be found in p4|, ^ and, within the numerical 
precision, it is given by 

< aiat(O) >= 1.000(1) . (5.5) 


On the other hand, the theoretical estimate of < (t( 0) > was given in eq. ( 3.11|) . Hence, 
comparing the two results, the relative normalization is expressed by the constant J\f as 


= M<j{x) = 0.930(3) a{x) . 


(6.6) 


Once these two quantities are hxed, there are no more adjustable parameters to com¬ 
pare the numerical data with the large-distance expansion of Gc{x). The form factors of 
the held cr(a;), entering the series (Q) can be easily recovered from those of 0(x), by 
using the relationship of these helds given by ( p.8|) , and for the correlation function we 
have ^ 

< a(a;)cr(0) >e= < 0(a;)0(O) >c • (5.7) 

The comparison between the two determinations of Gc{x) can be found in Figures 13 and 
14. In Fig. 13 we have only included the hrst three terms of Gc{x) (those relative to the 
form factors of the one-particle states Ai, A 2 , A^). As shown in this hgure, they can 
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reproduce correctly the behaviour of the correlation function on the whole infrared and 
crossover regions. A slight deviation of the theoretical curve from the numerical values is 
only observed for the hrst points of the ultraviolet region, where a better approximation 
can be obtained by including more terms in the form factor series. This is shown in Figure 
14, where hve more contributions (those relative to form factors up to state AiA^) have 
been added to the series. 


6 Conclusion 

The basic results of this paper can be summarised as follows. The Zamolodchikov S- 
matrix for the IMMF has been used as the starting point to implement a bootstrap 
program for the FF of the magnetization operator. Although the general solution of the 
bootstrap recursive equations remains a challenging mathematical problem, the matrix 
elements yielding the main contributions to the spectral representation of the correlator 
G{x) =< cr(a;)(T(0) > have been explicitly computed. This has enabled us to write a 
large-distance expansion for G{x) which is characterised by a very fast rate of convergence 
and provides accurate theoretical predictions for comparison with data coming from high 
precision numerical simulations. 

It would be interesting to obtain analogous results for the other relevant operator of 
the theory, namely the energy density e{x). To this aim, the only difficulty one has to 
face is the determination of the initial conditions for the form factor bootstrap equations 
appropriate for this operator. In the case of the held (j{x), we solved this problem by 
exploiting the proportionality with the trace Q{x) of the energy-momentum tensor. Notice 
that, due to the absence of symmetries in the space of states of the IMMF, the occurrence 
of the polynomials Qabi^^) the two-particle FF is precisely what is needed in order to 
distinguish between the matrix elements of o'(a;) and those of e{x). 

In conclusion, it must be remarked that the methods discussed in this paper can 
be generally used within the framework of integrable QFT. As a matter of fact, here 
they have been applied to a model which, for the absence of internal symmetries and 
the richness of its pole structure, can be considered as an extreme case of complexity. 
For instance, similar results to those contained in this paper can be obtained for other 
physically interesting situations, such as the thermal deformations of the tricritical Ising 
and three-state Potts models. The exact S-matrices for these models were determined in 


ref. 20, 35 . 
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Appendix A 


Aim of this appendix is to formulate in general terms the mathematical problem related 
to the computation of a generic form factor of the scalar operators in the IMMF. Such 
a formulation is obtained by exploiting a decisive property of the model, namely its 
bootstrap structure: any particle A, {i = 1 , 2 , ... 8 ) of the theory appears as a bound 
state of some scattering process involving the fundamental particle Ai and therefore can 
be obtained by a sufficient number of fusions of the particles Ai’s alone. The simplest 
examples are provided by the particles A 2 and A 3 , which appear in the initial amplitude 
S'ii( 6 '). Hence, the bootstrap structure of the model implies that all possible FF of the 
theory can be in principle obtained by the n-particle FF which only involve the particle 
Ai, by simply applying the residue equations (|3.5| ) the number of times we need to reach 
the FF under consideration. For instance, the form factors F22{0) and ^ 33 ( 6 *) can be both 
obtained by starting from Fim and by applying twice (|3.5|) on the poles at ^ and 
respectively. 

In view of the role played by the FF with the particles Ai, it is convenient to 
use a convenient notation. For brevity, we will denote them as F„( 6 *i, 6 * 2 ,..., 6 *„) = 
... ,9n)- It is now quite easy to hnd a parameterization of which correctly 
takes into account their monodromy properties and the pole structure. It can be written 
as 


AnjXi, ...,Xn) -pr 

{uJnixi ,... ,x„))” I}-. Dn{0ij){xi + Xj) 


(A.l) 


Let us explain the origin of each term entering the above equation. 

The monodromy equations ( |3.3| ) can be satished in terms of the functions 
solution of the equations 


nT"(9) = S„(9)F,T"(-9) . 

Fl7‘”(sir - 6 ) = J7’“”(sir + 9) . ^ ' 

These functions are required to have neither zeros or poles in the strip (0,27ri). 

can be explicitly written in terms of the functions G\{9) discussed in the appendix B as 

F;r(9) = -*smhjG|(9)Gj(9)G^(9) . (A.3) 

Once the monodromy properties of Fn are taken into account, we have to consider their 
pole structure. Note that, apart from the product of the the remaining part 

of these amplitudes can only be expressed in terms of functions of the variables 9ij which 
are even and 27ri periodic, i.e. functions of the variables cosh 6 'jj. Equivalently, they 
have to be symmetric functions of the variables Xi = e^\ A basis in the space of the 
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symmetric functions of n-variables is provided by the elementary symmetric polynomials 
uji{xi,X 2 , ■ ■ ■, Xn) 0 , defined by the generating function 

n n 

W{x + Xi)=^x'^~^UJj{Xi,...,Xn) . (A.4) 

k=l j=0 

The bound state poles of in all possible subchannels of the amplitude is encoded 
in the product of the terms 


Dn(0) = V.Mi)VMj)V.{e. 




(A.5) 


where V\{ 9 ) is dehned in eq. (^4.16|) . Concerning the kinematical poles, all of them are 
present in the product + ^i)- Finally, in ( A-l|) Hn is a normalization constant, 

A„(xi, X 2 ,..., Xn) is a symmetric polynomial and the last term (a;„(xi,..., x„))"' (which 
has no zeros in the physical strip) has been inserted in order to have a convenient form 
of the recursive equations. 

The polynomials A„(xi, X 2 ,..., Xn) can be obtained by solving the recursive equations 
(P^) and (|3.5|) . Using the parameterization ( |A.1|) , for the bound-state recursive equations 
we have 

A„+i(e“/^x, X2, • • •, x^) 


= A„(x,X 2 , ...,Xn) , 


(A.6) 


X ’^+3 [nfc=i(x + Xfc)(x - ?7Xfc)(x - 77 ^Xfc)] 
where r] = exp(77r/15). In writing ( |A- 6 | ) we have choose the normalization constants Hn 
as 


Hn+l — 


Tsin || sin ^ sin || sin 


2G'ii fcos? 


V 3 


COS T cos 


—V 

30 J 


sm 


2 llTT 
30 


47 


cos ^ cos -r COS 
3 5 


N ^ 

— I 
30 y J 


Hn 


(A.7) 


and, for simplicity, we have introduced the notation Gii{9) = G2(9) G2(9) G^{9). The 
constant 7 is given by 


7 = 


G2 (0)G2 I 

(¥) 

G2 

5 

(-f) 

IG2 1 

(f) 

IG^ 1 

15 

(-f) 

IG^ 1 

15 

(f) 

Gil 

(f) 

G2 

3 

(f) 

Gn (O)G^ (0)G2 (0)G_ 

-a(0) 

15 ^ ^ 


(A.8) 


A second hierarchy of recursive equations is obtained from the residue conditions of the 
kinematical poles, eq. (p.4|) , which can be written as 


{-)'^An+2{-X, X,Xi,...,Xn) = AnU (x, Xi, . . . , X„)A„(xi, . . . , X„) 


(A.9) 


where 


A = 


tg;, (^) (2 


COS o cos -r COS on 
3 5 30 


)' 


sin^ (r}i sin If sin ^ sin ff sin ff ^ 


2tt 


2tt 


7 sm — sm — sm 

8 Gii(0) sin"^ ^ 


(A.IO) 
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U{x,xi, 


)Xn) 2 ^ 


E (-) 


ki+k^+k^ 


X 


(A.ll) 


ki,...,kQ=0 


^6n-(fci+'"fc6) 


TT 


— (10(/ci — ^2) + 6(^3 — ^4) + (^5 — fcg)) 

L15 




Note that for spinless operators, the total degree of the polynomials A„ is equal to 
(3n^ — n)/2. The partial degree of A„ in each variable Xi is hxed by the asymptotic 
behaviour of the operator which is considered. For the FF of the spin operator a{x), the 
partial degree of A„ is given by (3n — 2). 

Given the initial values of the recursive equations, namely the one-particle FF Fi and 
the two-particle FF F2, one can solve eqs. ( A.ll ) and (|A-6|) in order to get higher-particle 
FF. For instance, for the three-particle case we have 


where 


Ag (3^4, X2, 3^3) — — ^\U}\U}2 A{u}\U} 2 F -\- Bu]-^U}2^'i 


2tt Ti 2 ti 

A = sm-h sm-h sm — , 

3 15 5 ’ 

27r TT 27r 27r TT 

B = 3 sin-h 2 sin-h 2 sin-h sin-sin — 

3 15 5 15 5 


(A.12) 


Appendix B 


In this appendix we discuss some useful formulas for the calculation of the functions 
pmin^Q'^ dehned in the text. Let us consider a function G\{9), solution of the equations 


Ga(0) = fxi9)G,i-e) ; 

Gx{i7r + 9) = Gx{i'K - 9) , 


where 


U9) = 


tanh ^ {9 + inX) 
tanh \{9 — inX) 


(B.l) 


(B.2) 


without zeros and poles in the strip (0,27ri). This function admits several different rep¬ 
resentations. One of them is given by 


Gx{9) = exp 


dt cosh 

- 1) \ 

2 9t 

Jo t cosh 1 sinh t 

27r 


(B.3) 


where 9 = in — 9. This function presents an inhnite number of poles and zeros outside 
the strip (0,27ri), as shown by the following inhnite-product representation 


cm = n 

fc =0 




1 + 


e/271 




1 , f 0/277 


1 + 


e/271 




fc+1 


(B.4) 
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or, equivalently 


gaw = n 

k=0 


r(| + A; + |)r(l + A;-|)r(l + A;+| + T (| + /c - I + 


r(| + A; + |)r(l + A; + |)r(l + A;-| + F (| + /c + | + 


(B.5) 

For numerical calculation, a quite useful expression is given by the mixed representation 


7V-1 

am = n 

A:=0 


1 + 


_e/2w 

fc+i-4 


1+(7%X 


1 I f g/2^ ^ ^ 

+ U+l+l 


1 + 


e/27i 


fc+l-l-l 


fc +1 


X 


exp 


foc dt cosh |(1 - 2A) 


/o t cosh I sinh t 


(iV + 1 - iVe-"‘)e 






From the above equations, we have 


and 


Gx{9) = G,.x{0) , 


Go{e) = Gi{e) = -isinh- 


The function Gx{6) satishes the following functional equations 

Gx{6 + i7i)Gx{9) = —i —^^-^(sinh^ + i sinTrA) 

sinTrA 


Gx{0 + i'K'y)Gx{9 - in'y) = 
G,.x{e)Gx-i{9) = 


'Gx{i7T'y)Gx{-i7T'jy 


Gx+^{9)Gx-,{9) 


Gx+-y{h)Gx-'y{0) 
sinh \ [9 — i(A — l) 7 r] sinh \ [9 + i(A + l) 7 r] 


sm 


2 ttA 
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(B.9) 
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Appendix C 

Aim of this Appendix is to illustrate through a specihc example the procedure used to 
obtain the results listed in Tables 3 and 4. Consider for instance the amplitude 822 ( 9 ). 
It exhibits hve simple positive residue poles located at 9 = 9 = 9 = 9 = 

and 9 = which correspond to the bound states Ai, A 2 , A 4 , A 5 and Ag respectively. In 
addition there is a pair of crossing-symmetric double poles at 9 = and 9 = The 
hrst one can be associated to the diagram of Figure 9.a with a = b = 2 and c = d = e = 1 
while the second one can be associated to the diagram of Figure 9.b with a = b = 2, 
c = d = f = 1 and e = 2 . A simple way to check that these diagrams actually satisfy the 
required constraints (all the particles on mass-shell and energy-momentum conservation) 
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is to draw the vector diagram for the momenta (dual diagram). Indeed, after performing 
the formal substitution 6 k —>■ il^k^ any two-momentum = {irtk cos /3fc, iruk sin j3k) can be 
thought as a vector of length ruk in the Euclidean plane. Then, the complete elasticity 
of the process requires that the external momenta form a closed parallelogram. Energy- 
momentum conservation also implies that the momenta entering each three-particle vertex 
form closed triangles; as a simple consequence of eq. (|2.8|) , the angle opposite to the side 
of length rric must be = in — If the triangles perfectly £t inside the external 
parallelogram, the original space-time diagram satisfy all the requirements. The dual 
diagrams associated to the double poles at 0 = and 6* = in the amplitude S 22 { 6 ) 
are shown in Figures 15.a and 15.b respectively. Of course, the internal parallelogram in 
the first diagram corresponds to the central interaction in Figure 9.a. 

It can be easily checked that 
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3z7r\ 

6 ~- ~ 

5 / 




U \2 


(r;,)"s„ (y) = (pj)'- 

According to the general discussion of Section 4, we write the FF F^{9) as 


(C.l) 

(C.2) 


F|(9) = (9) 




C22(9) 


where 


Notice that 


= -isinh- G 4 G 2 G 7 G 4 G 1 GI (6) , 

^■-4 > 2 L 5 3 15 T5 15 5J 


D22{6') = 


P4 P2 Pt_ P4_ P_i_ P3 P2 

5 3 15 15 15 5 5 


{ 0 ) , 


N 22 

Qg{9) = 4, cosh’^ 9 

k =0 



8 

T 

D22{9) ~ , 

i 
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(C.3) 

(C.4) 

(C.5) 

(C.6) 


In view of these asymptotic behaviours, eqs. (|4.9|) , (|4.10|) imply that N 22 < 3. Since 
F®, F^{9) and F®, a = 1,...,5 have been already determined in previous steps of 
the bootstrap procedure, eqs. (|4.21|) , ( [4.3D and ([4.31D provide the following seven linear 
equations for the four unknowns ®22; ^ — 0,..., 3 


F®(z7r) = 27iml 


F®{9 ~ zzz^a) - * 
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10 


lU 


c = 1,2,4,5 
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p2 t^2 

~ viliiil p0fol 


These equations turn out to be compatible and the solution is contained in Table 
particle Aq appears for the hrst time as a bound state in the amplitude 822 ( 0 )■ 
F^(9) has been fixed, can now be extracted using again eq. (|4.3|). 


(C.9) 

(C.IO) 

4. The 
Since 
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Table Caption 


Table 1 . S'-matrix elements of the IMMF where the numbers ( 7 ) in parenthesis should 
be read in units of They represent the functions tanh \{0 + / tanh \{0 — 

The bound state poles related to the particles Ai (i = 1,..., 8 ) are identihed by the 
numbers i placed above the functions ( 7 ). 


Table 2 . The hrst seventeen lowest energy states entering the spectral series. 

Table 3 . One-particle Form Factors of the operator 0. 

Table 4 . Coefficients of the polynomials Pah{0). 

Table 5 . The hrst eight contributions to the c-theorem. 

Table 6 . The hrst eight contributions to the universal amplitude of the free-energy. 

Table 7 . Numerical values of Gc{x) and their errors 6 for L = 64 and h = 0.075, as 
determined in [E3 . 





Figure Caption 


Figure 1 . Two-body S'-matrix. 

Figure 2 . Analytic structure of the S'-matrix in the s plane. The circles indicate the 
location of the bound state poles. 

Figure 3 . Bound state diagram in the s-channel. 

Figure 4 . Spectral representation series of the correlator. 

Figure 5 . Form Factor of the operator <h. 

Figure 6 . The two kinematical situations (a) and (b) entering the residue equations of 
the annihilation poles. 

Figure 7 . Bound state bootstrap equation for the Form Factor. 

Figure 8 . Residue equation for the bound state poles of the two-particle Form Factor. 

Figure 9 . Multi-scattering processes resulting in the double poles of the S'-matrix. 

Figure 10 . Residue equation for the simple poles of the Form Factor associated to the 
double poles of the S'-matrix. 

Figure 11 . Third-order pole in the S'-matrix. 

Figure 12 . Residue equation for the double poles of the Form Factor associated to the 
third-order poles of the S'-matrix. 

Figure 13 . Correlation function < (T(a;)(T(0) >c versus lattice space distances. The 
points on the graph represent the numerical data, as extracted from Table 7, while 
the continuum curve is the theoretical estimate with only the hrst three Form Fac¬ 
tors. 

Figure 14 . Correlation function < (T(a;)(T(0) >c versus lattice space distances. The 
points on the graph represent the numerical data, as extracted from Table 7 while 
the continuum curve is the theoretical estimate obtained with the hrst eight terms 
of the spectral series. 

Figure 15 . Dual diagrams of double pole graphs. 
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Ff 

0.4971505471 

Ff 

- 0.2627111760 

F3® 

0.1447685755 

Ff 

- 0.1107486745 


- 0.0467951944 

F® 

0.0336573286 

F® 

- 0.0127414814 

Fg® 

0.0023550931 
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1.623628945 
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